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Exposing the static scale of the glass transition by random pinning 
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Abstract. - The dramatic slowing down associated with the glass transition cannot be fully 
understood without an associated static length that is expected to increase rapidly as the tem- 
perature is reduced. The search for such a length was long and arduous, without a universally 
accepted candidate at hand. Recently a natural such length £ s was proposed, stemming from a 
cross-over between plastic and elastic mechanical responses of the material. In this Letter we show 
that supercooled liquids in which there exists random pinning sites of density p- lra ~ exhibit 
complete jamming of all dynamics. This is a direct demonstration that the proposed length scale 
is indeed the static length that was long sought-after. 
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A well accepted theory to explain the dramatic slow- 
ing down of the dynamics of super-cooled liquids upon 
approaching the glass transition is still missing, even af- 
ter decades of intense research efforts [HH]. One of the 
prominent approaches draws an analogy to the slowing 
down near continuous phase transitions; such an analogy 
requires a growing length scale that can accommodate the 
dynamical slowing down. The discovery of dynamic het- 
erogeneity in super-cooled liquids both in experiments [3J 
and in theoretical studies [2][5] and the detailed analy- 
sis of the associated dynamical length-scale using multi- 
point correlation functions [6ifl0] seemed to point in the 
right direction. But in Ref. [9] it was shown that the 
growth of the dynamical length-scale associated with the 
dynamic heterogeneity upon decreasing the temperature 
is not directly related to the rapid increase in relaxation 
time. It was also shown in Ref. [9] that the increase in 
relaxation time with decreasing temperature is better cor- 
related with the decrease in configurational entropy via 
the Adam-Gibbs relation [TT] which relates the structural 
relaxation time r a to the configurational entropy S c as, 
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Here A is a temperature independent constant. This phe- 
nomcnological relation is found to fit the data quite well 
over a large range of relaxation time in many experiments. 

The question of the existence of a static length that 
characterizes the degree of disorder in the approached 
glassy state remained open. Recently we were able to de- 



fine and measure such a static length-scale £ s [12 that 
grows upon approaching the glass transition. To be sure, 
there are other candidates in the recent literature, e.g a 
length-scale from the effects of boundary conditions, [15] . 
point to set correlations [14] , the scaling of the non affine 
displacement field [15] and patch correlation scale [16] . In 
this Letter we present strong arguments to demonstrate 
that the length-scale extracted by us via a different ap- 
proach [T2] is the relevant static length-scale in determin- 
ing the dynamics of the super-cooled liquids. 

The idea to single out the relevant length scale is based 
on introducing by hand another length scale, to observe 
the crossover phenomena expected in a system with two 
competing length-scales. If one has two length scales that 
affect a relaxation process then the dominant one will be 
always the smaller between the two. Thus by tuning the 
artificial length scale we can cross over from dynamics gov- 
erned by one length-scale to the other. A simple way to 
create a tunable length-scale is to introduce a density pi m 
of immobile paticles. The length-scale associated with this 
quenched disorder will be fj m ~ p^J A where d is the 
space dimension. If our proposed length-scale is indeed 
the static length scale then we expect to see a crossover 
from dynamics governed by £ s to a dynamics governed by 
the quenched disorder length-scale £j m . 

Before going into the details of the numerical analysis 
we should briefly mention how to extract the static length- 
scale £ s (details can be found in 12 ). Our starting point is 
the fact that at low frequency tail of the density of state 
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Fig. 1: (Color online) Left panel : The measured minimal eigenvalue as a function of system size for a binary system interacting 
via Kob- Andersen potential in 3 dimension. Every data point represents an average over 2000 Inherent structures which are 
the local energy minimum of the potential energy function. Right panel : Same data set but plotted according to the scaling 
ansatz Eq. [3] to obtain the static length-scale. Inset shows the temperature dependence of the obtained length-scale. 



(DOS) of amorphous solids consisting of N particles re- 
flects the excess of plastic modes which do not exist in the 
density of states of purely elastic solid. [TTlfTB] . This excess 
of modes is sometime referred to as the 'Boson Peak' [15] . 
Here and below the 'mode' refers to the eigenfunction of 
the underlying Hessian matrix. Recently |20j we discov- 
ered that the eigenvalues {Xi}i — l dN (with d being the 
space dimension) appear in two distinct families in generic 
amorphous solids, one corresponding to eigenvalues of the 
hessian matrix that are only weakly sensitive to external 
strains; the other of eigenvalues go to zero at certain val- 
ues of the external strain, thus leading to a plastic failure. 
The first group of modes is decently described by the the 
Debye model of an elastic body, but this is not the case for 
the second group corresponding to the density of plastic 
modes. 

For the purposes of the present letter it is enough to 
write this excess part (in the thermodynamic limit) as 
B(T)j pi (-^y) > where the pre- factor B(T) being strongly 
dependent on temperature. Particular models for this part 
were presented in [31]. We do not need to specify a func- 
tion here, and it is only important to understand that 
this function is a partial characterization of the degree of 
disorder which grows upon approaching the glass transi- 
tion. Together with the standard Debye contribution one 
can approximate the low-frequency tail of the density of 
states as 
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Here A^ ~ /xp 2 ' d_1 , is the Debye cutoff frequency and \jl 
is the shear modulus. The physical idea that allows the 
determination of the static typical scale is that the min- 
imal eigenvalue A m ; n observed in a system of N particles 



will be determined by either the first or the second term in 
Eq. ([5]). For a system large enough, local disorder will be 
irrelevant in determining A m j n , and it will be decided by 
the Debye contribution. For small systems the opposite is 
true. Thus there exists a value of N where a cross-over 
occurs. This cross-over is interpreted in terms of a typical 
length-scale separating correlated disorder from asymp- 
totic elasticity. 

In Ref. [12j it was shown how to derive an implicit equa- 
tion for where the angular brackets represent an en- 
semble average over many systems prepared at the same 
temperature. The equation takes the form 



(A m in) 



T 



Hi?) 




(3) 



The typical scale £ S (T) will be calculated by demanding 
that all the data calculated for different system sizes and 
temperatures should collapse into a master curve just by 
appropriately choosing the £ S (T). 

As an example consider the analysis done for the glass 
forming Kob- Andersen binary mixture [22] at number den- 
sity p = 1.20. The systems were equilibrated at some 
temperature T > T g and then inherent structures are cal- 
culated by direct energy minimization to the nearest lo- 
cal minimum of the potential energy landscape. At this 
state the Hessian was computed and the minimal eigen- 
value was obtained using the Lanczos algorithm [23]. For 
a given system size N and temperature T this procedure 
was repeated to have an average (Amin) until convergence 
was achieved. At this point the temperature or the sys- 
tem size were changed and the procedure was repeated, to 
eventually have a table of (X m i n )(N,T). In the left Panel 
of Fig. [U we have plotted the minimal eigenvalue rescaled 
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by the characteristic Debye value as a function of system 
size N for different temperatures. In the right Panel the 
same data is plotted according to the scaling ansatz Eq. Q 
to extract the length-scale £ s by collapsing the data into 
a master curve. The resulting data collapse into a master 
curve indicates that our scaling ansatz Eq. © is obeyed 
to high precision. In the inset at the right panel we show 
how the typical scale increases when the glass transition 
is approached. 

We also showed in [12] how our typical scale helps in 
understanding other measures of disorder that were pro- 
posed in the past. As an example we showed that the 
system size dependence of configurational entropy S C (T) 
[9 , 24; (the reader is referred to these publications for a 
full description of the method and the results) of the Kob- 
Andersen model in 3 dimensions can be explained using 
the typical scale proposed in [12] . This also gives us a 
direct relation between structural relaxation time and the 
static length-scale using the Adam-Gibbs relation [11] . 
Below we will see how this relation along with the static 
length-scale helps us to understand the effect of immobile 
particles on the dynamics of supercooled liquid. 

We start with the Adam-Gibbs relation Eq. (JTJ. In the 
presence of immobile impurities of density pi m the config- 
urational entropy changes, and we propose that it has a 
scaling form, reading 



S c (p im ,T) = S c (T)g(p im ^(T)) 



(4) 



where g(x) is an unknown scaling function. Thus with the 
presence of randomly pinned particles the relaxation time 
will depend on the scaling function g(x) and Pi m S,f as 



T a (pim,T) = T exp 
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The scaling function g{x) must have the following asymp- 
totic behavior 



g(x) — > 1, as x — > 0. 



(0) 



So using Eq. ([5]) and Eq. ([T]), we end up with following 
equality 



log 



7~a (_Pim ? -^J 


= A(1 


. t q (0,T) _ 





- 1 



(7) 



Now using Eq. we can expand the scaling function 
g{x) for small x as 



g(x)~l-\g \x + O{a?), 



(8) 



We expect g(x) to be a decreasing function of x as we no- 
ticed before that with increasing disorder density pi m the 
relaxation time T a increases so the first derivative calcu- 
lated at x = 0, \g \ will be negative. Putting Eq. (|S]) in 
Eq. Q we arrive at the following relation after little bit 
of algebraic manipulation, 
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Fig. 2: Average Overlap function (Q(t)) is plotted for different 
disorder density p im for T = 0.700. 



This relation is expected to be valid only in the dilute dis- 
order limit. Next we explain the details of the experiment 
and the results that validate this scaling theory. 

The experiment is done as follows. We first equilibrate 
a system of N = 1000 particles interacting via the Kob- 
Andcrsen Potential [22] . The interaction potential is given 
by 



,(T. 



V aP (r)=4e a0 [C-^) 12 ~C-^) e } 
r r 
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where a,0£ {A, B} and = 1.0, e^s = 1.5, ess = 0.5, 
<j aa = 1.0, <tab = 0.80, (tbb = 0.88. The Interac- 
tion Potential was cut off at 2.50a a p. We have per- 
formed the simulations at six different temperatures in 
range T e [1.00, 0.45]. Then we take the equilibrated 
system and randomly choose N im — N * p irn number of 
particles and freeze their degrees of freedom and run the 
dynamics to calculate the relaxation time. 

We have calculated the overlap function defined below 
to estimate the structural relaxation time r a , 
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where the weight function w(x) — 1 if a; < 0.30 and zero 
otherwise. The prime sign in the summation indicates that 
immobile pinned particles are not included in the summa- 
tion while calculating the correlation function. In Fig. [21 
we plotted the overlap function Q(t) averaged over 50 dif- 
ferent realization of quenched disorder for different density 
of disorder at temperature T = 0.700. The Relaxation 
time t q is defined to be the time where the correlation 
function Q(t) — 1/e, where e is the base of natural loga- 
rithm. 

In Fig. [3] we have shown the dependence of T a as a 
function of disorder density pi m for the six simulated tem- 
peratures. Note that as one decreases the temperature the 
effect of disorder in the relaxation time kicks in for smaller 
disorder density. This can be understood easily from the 
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Fig. 3: The relaxation time r Q measured for all temperatures 
plotted against density of pinned particles pi m - 



fact that the effect of the quenched disorder is to pin the 
correlated region to one place causing hindrance to relax- 
ation. Now as one decrease the temperature the size of 
these correlated region grows and one needs less number 
of pinning sites to freeze the whole system. 

In Fig. 21 we have plotted TS C log [r a (pi m ) / V a (0)] mea- 
sured for all the simulated temperatures against pi m and 
in Fig. El we have plotted the same data but rescaling the 
abscissia by the appropriate power of static length-scale £ s . 
The quality of the collapse suggests that indeed our pro- 
posed scaling ansatz Eq. (J7]) is obeyed to a high precission. 
The solid line is the approximate scaling form derived in 
Eq. ((9]). It is clear from the figure that the approximate 
form of the scaling function is indeed a very good approx- 
imation in the dilute disorder regime. We also want to 
point out here that Eq. © suggests that one can have an 
ideal glass transition as a function of disorder density at 
a critical density of disorder p\ m (T) given by 



Pim(T) 
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(12) 



In summary, we have showed that the dynamics of su- 
percooled liquid under the influence of externally imposed 
quenched disorder can be understood completely by our 
proposed static length-scale. This analysis also suggests 
that one can get the static length-scale by studying the 
dynamics of the supercooled liquids under externally im- 
posed quenched disorder instead of calculating it from the 
scaling of minimal eigenvalue which may be less trivial to 
access in experiments. It would be really nice if one could 
experimentally study the dynamics of supercooled liquids 
under externally imposed disorder, since the configura- 
tional entropy which is needed to have the data collapse 
is also easily calculable in experiments. We hope that this 
kind of analysis will be done on other model glass form- 
ers and also in laboratory experiments by other groups to 
enhance the understanding of the glass transition. 

This work had been supported in part by an ERC 



Fig. 4: The logarithm of the rescaled relaxation time 
T a (pim)/Ta(0) multiplied by TS C measured for all the simu- 
lated temperatures plotted against p.; m . 
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Fig. 5: The data collapse of relaxation time using the scaling 
ansatz[7] The solid line is the approximate scaling form derived 
in Eq. M 
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